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Abstract. We prove that the set of solutions to the paraboUc singular p- 
Laplace equation with Dirichlet boundary conditions on a bounded Lipschitz 
domain O for all space dimensions is continuous in the parameter p £ [1, +00) 
and the initial data. The highly singular limit case p = 1 is included. In 
^>f^ , particular, we show that the solutions Up converge strongly in (Q), uniformly 

in time, to the solution ui of the parabolic 1-Laplace equation as p — > 1. 



O 



', 1. Introduction 

^ ■ Let n C W\ d € ¥1 he an open bounded domain with Lipschitz boundary 

I dfl. Consider the following family of singular diffusion equations indexed by p G 

^' (l,+oo), 

^ . f ^ - div [\Vup\P-^Vup] = fp in (0, +00) X n, 

> '. (pd) S Up = on (0, +c5o) x dQ, 

^ . [ Up{0) = Xp, 

C<i '. where t.^ P ^?(Q^ arid P T,? CCn. 1,^ 



O 



■ Consider also 

O 



where Xp S L^(ri) and fp £ Lf^^{{0, +oo)-, L^{^)). The notation (po) refers to 
Dirichlet boundary conditions. 



^-div[Sgn(VMi)] 9/i in (0,+oo)xr!, 
Ml = on (0, +00) X dfl, 

Ml(0) = Xl, 



^ . where xi G L^(f2) and /i G Lf^^{{0, +00); L'^{il)) and the vector-valued (and multi- 

. .' valued) sign-operator Sgn : IR"* ^ 2^ is defined by 



Sgn(a;) 



if a; G R'' \ {0}, 



\x\ 

{2/ e R"* I |y| < 1} , ifa; = 0. 

In fact, the expression "div [Sgn(V-)]" is merely heuristic, we shall later replace it 
by the subdifferential c^'I'i, where the convex potential $1 is defined in (2.6) below. 

In this work, we consider weak variational solutions to equation (pd), P G 
[l,-|-oo), defined as limits of Yosida-regularized equations. We are interested in 
continuity properties of (pd) in the initial data and, in particular, when the pa- 
rameter p varies. The p-Laplace evolution (pn) belongs to the class of gradient 
flow-type equations, where the dynamics is generated by the (infinite dimensional) 
gradient of a differentiable functional. 
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In order to prove the desired continuity of solutions in p, we shall apply a the- 
orem by Attouch, stating that so-called Mosco continuity of convex, lower semi- 
continuous potentials implies the continuous dependence of the solutions to the 
related, subgradient-type evolution equations, see [6]. The classical theorem re- 
quires a Hilbert space framework, which explains our choice of L^-data. For Banach 
spaces, variational convergence and viscosity solutions are brought together in [8]. 

Inclusion (Id) requires some attention. It is called total variation flow, see [3,4]. 
We shall model it by the gradient of the total-variation functional on functions of 
bounded variation in L^(r2). 

We shall also consider the p-Laplace evolution with Neumann boundary condi- 
tions, 



(pn) 



011 

- div [I Vup|P-2VuJ = fp in (0, +00) X n, 



dup 



dv 

Up{0) = Xp, 

where v is the generalized outward normal on dO.. 
Also 



on (0, -l-oo) X dVt, 



(1^) 



^-div[Sgn(Vni(i))] 9/i(t) in (0,+oo)xO, 

=0 on (0, -l-oo) X dVl, 

ov 



■"1(0) 

where again the expression "div [Sgn (V-)]" is made rigorous as a subdifferential 
d^i of the potential defined in (2.7). 

The singular diffusion operator in equations (pd): {pn) is called p-Laplacian. 
Both equations of evolution-type are covered by the theory of nonlinear semigroups 
associated to equations of subdiff'erential-type in Hilbert space as in the book of 
Brczis [17]. For the notion of solutions used in this paper, sec Definition 3.2 below. 
Explicit representations for the subdifi'erential operators were given by Schuricht 
[39] for inclusion (l^) and by Andreu, Ballester, Caselles and Mazon [3] for inclusion 
(Ia^)- 

Classically, the p-Laplace equation appears in geometry, quasi-regular mappings, 
fluid dynamics and plasma physics, see [21,22]. In [28], Ladyzenskaja suggests the 
p-Laplace evolution as a model of motion of non-Newtonian fluids. A typical 2- 
dimensional application can be found in image restoration, see [13, Ch. 3] for a 
comprehensive treatment. General equations of p-type are studied in [41]. 

Eigenvalue problems for the 1-Laplacian have been studied e.g. by Pridman, 
Kawohl, Schuricht and Parini [25,27,37,39]. 

The stochastic analog to (po) is being examined by Liu [32], respectively, to (Id) 
by Barbu, Da Prato and Rockner [15]. Recently, in the case of stochastic equations, 
similar results to this paper have been proved by Ciotir and the author in [19]. 
The proof in this paper, however, originates from the author's thesis [40, Ch. 8.3], 
where variable-spacc-approach, similar to that in [11], is used. Previous convergence 
results for the stationary problem have been obtained by Mercaldo, Rossi, Segura 
de Leon and Trombetti in [34-36] . Convergence of the evolution problem for local 
solutions has been previously investigated by Giga, Kashima and Yamazaki in [26] . 

Let us formulate the main result of this paper. Its proof can be found in section 
4 below. The notion of solution is specifled in section 3 below, see Deflnition 3.2. 

Theorem 1.1. Let po € [l,-|-oo), {pn} C (l,+oo) with lim„p„ = po- Let T > 0. 
Suppose that lim„ \\xp^ - Xpjj^^^^^^ = and lim„ ||/p„(i) - /po(*)IIi,2(q) dt = 0. 
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Denote by Up the solution to equation (pa) in the sense of Definition 3.2 and by Vp 
the solution to equation (pn) Definition 3.2. 
Then 

(1.1) lim sup -WpoWllL^m) =0 
and 

(1.2) lim sup \\Vp„{t)-Vp„{t)\\^^2m)=0- 

" te[o,T] ^ ' 

In other words, if 

F : [1, +oo) X L\n) X Ll^{{0, +oo); L^fl)) ^ C{[0, +oo); L^fl)) 

X f 

denotes the map that assigns the solution 1 1-> Up"' "{t) to equation {pd), our result 
states that F is continuous. Certainly, a similar statement holds for (pn)- 

2. The related energies 

We observe that, for p € {1, +oo), the p-Laplace operator on L^(f2) is a gradient- 
type operator with the convex potential 



%{u) := 



i^j |Vu|Pdx, ifi 
l^+cx), ifi 



(2.1) 



ifue Wo'P(0)ni2(0), 
^uei2(0)\Wo'f(n). 

Note that the Dirichlet boundary conditions are encoded into the energy space 
Wo^'^(f2) which is the standard first order p-Sobolev space on fl such that its ele- 
ments have vanishing trace (see next paragraph). 

In general, for an element u G W^'P{CI), we define its trace ^p{u) to dO. by 

Ion 

= I Mdivr;da;+ / {r],Vu) dx \fr] & C^{n;]R'^), 

where ly is the generalized outward normal on dil and is the d— 1-dimensional 

Hausdorff measure on d^. We have that 7p : VFi'P(f^) LP{dn,M"^-^) is a 
continuous linear operator. Note that our definition of Wq'^{^) coincides with the 
standard one, where one takes W^'^iVt) to be the closure of C^{Vl) in W'^'P{Vt), 
see [23, §5.5, Theorem 2]. We denote by di, V resp. the weak or ordinary partial 
derivative in direction of the i-th coordinate, the weak or ordinary gradient resp. 
Di and D refer to the corresponding objects in the sense of Schwartz distributions. 
Sometimes we shall write [Df] instead of Df to indicate that [Df] is a Radon 
measure. Of course, if / e L\^X^) ^"^^ ^ -^ioc(^; ^^^n D/ = V/ a.e. 

Define also the energy of the p-Laplace operator in L^(0) with Neumann bound- 
ary conditions by 

.\pdx, iiuew^'P{n)nL^{n), 
Hue L^{n)\w^'P{n). 

Let us continue investigating It is easily seen that $p is a proper convex 
function. It is well-known that $p is lower semi-continuous, see e.g. [18]. 

The corresponding subdifferential d^p is a realization of the p-Laplace operator 
on Q with Dirichlet boundary conditions. Recall that, in general, for a convex, 
proper, lower semi-continuous function ^ : H ^ (— oo,-|-oo] on a separable Hilbert 
space H, the subdifferential d^, defined by [x, y] G 9^* iff 

{y, z-x)h< ^{z) - *(a;) e H, 



III 

[+0O, 



4 



J. M. TOLLE 



is a maximal monotone graph. 

On smooth functions G Cq°(0), has the representation 

{d%){^) = -div [|V<^|f-i sgn(V<^)] , 

where sgn is the selection of Sgn with minimal Euclidean norm. 

Note that the extension of <I>p to W^'^ifl) is Frcchet differcntiable in ^^'^(il)- 
norm, a fact which we, however, will not use. By the theory of subdifferential 
operators, dom(9$p) is dense in 1/^(0), see e.g. [14, Ch. 2, Proposition 2.6]. 

We shall continue investigating the limit case p = 1. We need some facts from 
the space BV, which we shall collect here. 

Definition 2.1. A Junction f e Ll^^{il,) is said to be of bounded variation if 
IID/II (O) :=sup| / /divr?da; 



The value ||D/|| {ft) is called the total variation of f. The space of all classes of 
functions in L^{Sl) that are of bounded variation is denoted by BV{fl). 

Let / £ BV{rt). Then there is a Radon measure pi f on fl and a measurable 
function af : Q ^ Wi'^ such that |cr/| = 1 /U/-a.e. and 

(2.2) f fdiyr]dx = - f {r,,af) dpif Vt? e Co^(n; R'^). 
Jn Jn 

By [24, §5.1], ||D/|| (O) = /i/(0). Set d[D/] := afiJ,f{dx), which is a IR'^-valued 
Radon measure on O. Hence (2.2) becomes 

(2.3) / /divr?dx = - / (r?,d[D/]) Vr? e Cq^O; Il<^). 

For each / e BV{il) there is the trace ji{f) to dil satisfying the following formula 

(2.4) / /divrydx = - / (7?,d[D/])+ / j,{f) {r],iy) dJ^"-' Vry e ^1(11; R"^), 
Jn Jn Jdn 

where v is the outward normal and ^f^^^ is the Hausdorff measure on d^. We 
have that 7i(/) e L^{d^,M"^-^). We refer to [24, Ch. 5.3] and [2, Ch. 3] for 
details. 

Also, we have that ^"^'^(0) C BF(0), p G [1, +oo] and 

(2.5) /|Vu|da;= ||Du||(f]) Vu G VF^'P(r2), 
Jn 

see [24, Ch. 5.1, Example 1]. 

Remeirk 2.2. Suppose that dil is Lipschitz. Let u G BV{fl) and extend u by zero 
outside O. Then u G BV{W'-) and 

||D«|| (R^) = ||Du|| (17) + / \^r{u)\dje''-' 

Jdn 

where 71 (u) G V-{dQ.,M"'--'^) is the trace of u. We refer to [2, Theorem 3.87]. 

We are ready to define the convex potential associated to the 1-Laplacian in 
L^(il) with Dirichlct boundary conditions. The Dirichlet problem was derived by 
Kawohl and Schuricht in [27,39]. 

Define 



(2.6) $i(u) := 



||Du|| (R<^), if u G BV{Q.) n L2(0), 
+00, if uGi2(n)\W(0). 
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$1 is obviously proper and convex. Lower semi-continuity follows by standard 
arguments via integration by parts since ||Du|| (R'') is lower semi-continuous with 
respect to weak convergence in Ljq^(R''). 

The 1-Laplacian is defined to be the subdifferential 9<i?i of As above, 

dom(9$i) is dense in L^{fl). A precise representation of is given in [39, The- 
orem 2.1] and [27, Theorem 4.23]. 

Define also 



T)u\\{n), a ue BV{n)nL^{n), 

-00, if w e L'^{n)\BV{Q). 



(2.7) *i(u) := 

\E'i is convex, proper and lower semi-continuous, see [4] 



3. MOSCO CONVERGENCE 

Let H he a separable Hilbert space. Let : — > (— oo, +oo], n e N, $ : if — )• 

(— oc,+oo] be proper, convex, l.s.c. functionals. 
Recall the following definition. 

Definition 3.1 (Mosco convergence). We say that $ in the Mosco sense if 

(Ml) \fx e H \fxn G H, n e'N, Xn ^ X weakly in H : liminf $„(a;„) > $(x). 



(M2) \/y e H 3yn € H, n e N, y„ y strongly in H : limsup$„(y„) < <^{y). 

n 

Consider the following sequence of abstract gradient-type evolution equations in 

H: 

^3 I ^Un{t) + d^n{Un{t)) 9 /„(t), 0<t< +00, 

[ U„(0) = Xn, 

where Xn & H and /„ e L'^^^{{Q,+cc);H). 
Also, 

(3 2) {^u{t)+d^{u{t))^f{t), 0<t<+oo, 

1 «(0) = X, 

where x € H and / e iu)c((0' +0°); H). 

A solution to, say, equation (3.2) is defined as follows. Compare with [17, 
Definition 3.1]. 

Definition 3.2. We call u G C([0, +oo); if) a solution to equation (3.2) if u is lo- 
cally absolutely continuous, u(0) = x, u{t) e dom(9$) for a.e. t > and (3.2) holds 
for a.e. t > or, equivalently, in the sense of distributions in Lf^^{{0,-\-oo);H). 

We shall need the following theorem. 

Theorem 3.3. Suppose that for T > we have that a;„ — )• a; strongly in H, fn ^ f 

strongly m L'^{{0,T); H). 

Then $„ <& implies that Un — > u strongly in H, uniformly on [0,T]. 
Proof See [7, Theorems 3.66 and 3.74]. □ 
See [5,6,9,20] for related results. 
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4. Proof of Theorem 1.1 

Let po G [l,+oo), pn G (l,+oo), n € N such that hm„p„ = po- Write g„ := 
Pn/{Pn — 1), n G N, qo := po/(po — 1) (with 1/0 := +oo). If we can prove that 

(4.1) ■■= ^ > *po =: 

n— >+cx) 

the assertion of Theorem 1.1 follows from Theorem 3.3. 

Before we prove (4.1), we shall need the following approximation result. 

Lemma 4.1 (Littig-Parini-Schuricht). Let O C M.'^ be open and bounded with 
Lipschitz boundary dil and let u € BV{^) fl L^'(n) for some p € [l,+oo). Then 
there is a sequence {uk} in Cg^{Q) such that, for any p G [l,p], 

(4.2) \\uk - ^ and \\Buk\\ (R'^) ^ ||D7i|| (R''). 

The result was proved by Littig and Schuricht in [31, Theorem 3.2], see also [30, 
Satz 2.31] . Previously, if dCl is C^, the result was proved by Parini, see [37, Theorem 
2.11] and also [38]. 

Proof of (4.1). We start with proving (Ml) from Definition 3.1. Suppose that 
Un e -L^(ri) with 

liminf $„(w„) < +oo. 

n 

Extract a subsequence (also denoted by {u„}) such that 

lim$„(u„) = liminf $„(u„) 

n n 

and 

C := sup$„(u„) < +00. 

n 

Since m„ G Wq'^" (fl) and hence '^p„{un) = 0, we can extend u„ to R'^ by zero 
outside fl (denoted also by u„) and get that u„ G W^'^"(R'^) (cf. [29, Exercise 
15.26]). 

Fix a ball B c R''. Suppose first that for a subsequence we have that < 2. 
Then by Holder inequality. 




and 

/ \WUn\dx < (P„$„(W„))'/^" (volS)l/«" 

(4.3) Jb 

< (l + Csupp„)(l+volB) 

n 

Suppose for a while that 

(4.4) sup ||wn||z,2(]a<i) < +00. 

n 

Hence is bounded in W^'^{B). If for a subsequence Pn > 2, the W^'^-bound 
can be established with the help of Jensen's inequality. 

By compactness of the embedding W^'^{B) C L^{B), a subsequence of 
converges in L^{B) (see e.g. [33, §1.4.6, Lemma] or [1, Theorem 6.2]). By a diagonal 
argument, we can extract a subsequence such that converges strongly to some 
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/ e Ll^^CR''-). W.l.o.g. Un ^ f dx-a.e., by extracting another subsequence, if 
necessary. Also, the measures {diUn dx} are vaguely bounded and hence vaguely 
relatively compact, see [16, Paragraphs 46.1, 46.2]. For each 1 < i < d, we can 
extract a subsequence, such that {diUn dx} converges to some locally finite Radon 
measure rrii on R''. By vague convergence and integration by parts, 

y i^dmi = lim J cpdiUn dx = — lim J di(pundx = — J dipfdx, 

for every if € C§°{]R.'^) and every 1 < i < d. Hence mj = Dj/. Furthermore, for 
every if G C^(R''; R''), 



/ 



Un div ip dx 



1 

Pn 



(4.5) 



<$„(m„) X < 



f ||<^|U(vol(supp<^))^"/* 

Pn/Qn 

da:' 



{VUn,p) dx 

if Po = 1, 
if Po > 1- 



Upon taking the limit, by Lebesgue's dominated convergence theorem (since either 
for a subsequence |-|^" < 1q + |-p or for a subsequence |-|'" < la + |-|^"Pn9") 



(4.6) 



— f f div if dx 
Po J 

< liminf $„(««) x < 



[(/' 



)Po/?o 



if Po = 1, 
if Po > 1- 



Taking the supremum over all if with \\ip\\^ < 1 (if Po = 1) or with ||<p||i9o(]R<i) ^ 1 
(if Po > 1) yields 

IID/II (R<^) if PO = 1, 
1 



Po 



ID/If" da; ifpo > 1 



< liminf i>„(u„). 



Suppose now that u„ ^ u weakly in L'^{{1). This justifies (4.4). Clearly, for all 

ip e Co(R'^), 

J uipdx = lim J Un<pdx = J fipdx, 

hence u = f dx-a.e. and Du = D/ in the sense of distributions. We are left to 
prove that, if po = 1, then u € BV{ft), because then 

llDull (R'') = $(w) < +00. 

Simiarly, if po > 1 and if m £ Wo'P°{n), 



$(u) = — / \Vu\P° dx< — [ iDuf dx < +00. 
Po Jn Po J^d 



We have that u„ — > u in L\^^{W^) and that diUndx DjU in the vague sense 
on R''. Hence by the definition of the trace (2.1), (2.4), 

= lim / 7p„ [un) {ip, y) dM"^-^ = [ 7j,, (u) {ip, u) dM'^-^ 
" JdQ. Jan 

V(^G Co°°(R'';R'^). 

Let ip e C^{n;W^). By (2.4), 

/ udiviy9da; = — / (lyj, d[Du]). 
Jn Jn 
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Hence by definition, u £ BV{^). 

On the other hand, if po > 1, we can take the the supremum over all G 
C^{Wi'^;W^) such that ||¥'||i,.o(a) < 1 to get that 

/ |7po(«)r° dJf"^"' < liminf / IIp^K)^" ^^"^'^ = 0. 
Jon " Jdn 

Compare with [2, Lemma 3.90] and (4.6). We get that 7po('«) = Jf"^~^-a.e. and 
hence u € Wq'^°{Q). Since we can repeat the steps for any subsequence of {un\, 
we have proved (Ml). 

Let us prove (M2) from Definition 3.1. Let pq = 1 and u G BV{Vl). Then by 
Lemma 4.1 there is a sequence {um} C Co°(ri) C BV{Q,) with 

u„ ^ u in L'^{n) and ||Du„|| (R'') ||Du|| (R'') as m ^ +oo. 

But by Lebesgue's dominated convergence theorem (since for large n it holds 
that ^H^" < 1q + \f) and (2.5) for each m e N 



— / |Vum|^" da; -)■ / |Vuto| da; = ||Duto|| (R**) as n -)■ +oo. 

An application of Lemma A. I shows that there exists a sequence {nin}-, rUn t +oo, 
such that 

lim$„(«„J = lim— / |VM„Jf" da; = ||Dw|| (R'') = $(«), 

which proves (M2). If po > 1> we can repeat the above steps in the obvious 
manner. □ 

Let us briefly treat the case of Neumann boundary conditions. The idea of proof 
is essentially the same with slightly different methods. Therefore, let po G [1, +oo)) 
Pn G (l,+oo), n e N such that lim„p„ = po- Write := Pn/{Pn — 1), n G N, 
90 := Po/{po — 1) (with 1/0 := +oo). In order to apply Theorem 3.3, we would like 
to prove that 

(4.7) := =: 

We need the following approximation result. 

Lemma 4.2. Let u G BV{Q.) fl L^(0). Then there exists a sequence {uk} C 
C°°(n) n BV{n) such that 

11'"'= -"IIl^co) o-'^d ||Dufe|| (Jl) ||Du|| (O). 

Proof. With replaced by this is well-known, sec e.g. [10, Theorem 10.1.2]. 
Since the approximation is obtained by a mollifier, we get the strong convergence 
in Z/^ by Hausdorff- Young inequality. □ 

Proof of (4.7). We start with noting that by the Lipschitz assumption for dVL, and 
the fact that W^''P{n) C BF(fi), p G (1, +oo), wo can sot any function u G W^^'P{n) 
a.e. zero outside (denoted also by u) and obtain a function u G BV{Uf'), see [24, 
§5.4, Theorem 1]. 

Let us prove (Ml) from Definition 3.1. Suppose that u„ G Z/^(0) with 

liminf ^'„(u„) < +oo. 

n 

Suppose also that ii„ u weakly in L^(0) for some u. Extract a subsequence (also 
denoted by {wn}) such that 

lim*„(u„) = liminf *„(u„) 
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and 

C := sup < +00. 

n 

Hence u„ G W^'P"{Q). Extend u„ by zero outside Cl, denoted also by m„. Note 
that M„ e L'^QR!^) and w„ e ST/(]R'*), but possibly u„ ^ ^^■^"(R'^). 

As above we can extract a subsequence of and some / e BV{il) such that 

^ / in Lj'Q^(ri) and the distributional gradient of converges to some locally 
finite Radon measure m which is then found to be equal to [Df]. In fact, Un ^ f 
in L'^{n) by [24, Theorem 5.2.4]. 

The rest of the proof works exactly as above except that the test-functions in 
C^(]R'*; R'^) are replaced by test-functions in (fi; R'') and in (4.6) we take the 
supremum over all (p such that ||'^||j;,5o(n) ^ 1 instead. 

The proof of (M2) from Definition 3.1 can easily be completed by the arguments 
of the above proof combined with Lemma 4.2 resp. the well-known approximation 
of W^'P {nyfunctions by C°°{n) n 'J' (^^)- functions, see e.g. [24, §4.2, Theorem 
3]. □ 
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Appendix A. A diagonal lemma 

Lemma A.l. Let {a„,m}ra,meN C R 6e a doubly indexed sequence of extended real 
numbers. Then there exists a map n i->- m(n) with m{n) t -l-oo as n ^ +00 such 
that 

(A.l) liminfa„TO(n) ^ liminf 

n— >-t-oo ' m— f+cx) 

or, equivalently 

(A. 2) lim sup a„_m(n) < limsup 

Proof. See [12, Appendix] or [7, Lemma 1.15 et sqq.]. □ 



lim inf o„ 

n— >-t-oo ' 



lim sup an,n 

n— ^+00 
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